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ABSTRACT 

We observed the gamma-ray pulsar Geminga with the FUV-MAMA and NUV-MAMA detectors of the Space 
Telescope Imaging Spectrometer to measure the Geminga's spectrum and pulsations in the ultraviolet. The 
slope of the far-ultraviolet (FUV) spectrum is close to that of a Rayleigh- Jeans spectrum, suggesting that the 
FUV radiation is dominated by thermal emission from the neutron star (NS) surface. The measured FUV flux, 
Fpuv = (3.7 ± 0.2) x 10~ 15 ergs cm" 2 s" 1 in 1 155-1702 A band, corresponds to a brightness temperature 7rj w 
(0.3-0.4) (d2oo/R\3) 2 MK, depending on the interstellar extinction (d = 200^200 pc and R = 13/?b km are the 
distance and the NS radius). The soft thermal component of the Geminga's X-ray spectrum measured with the 
XMM-Newton observatory corresponds to a temperature T s = 0.49 ±0.01 MK and radius R s = (12.9 ± 1.0)c?200 
km. Contrary to other NSs detected in the UV-optical, for which the extrapolation of X-ray thermal component 
into the optical underpredicts the observed flux of thermal radiation, the FUV spectrum of Geminga lies slightly 
below the extrapolation of the soft thermal component, which might be associated with the Geminga's very low 
temperature. Surprisingly, the thermal FUV radiation is strongly pulsed, showing a narrow dip at a phase close 
to that of a broader minimum of the soft X-ray light curve. The strong pulsations might be attributed to partial 
occultations of the thermal UV radiation by regions of the magnetosphere filled with electron/positron plasma. 
In contrast with the FUV spectrum, the near-infrared (NIR) through near-ultraviolet (NUV) spectrum of 
Geminga is clearly nonthermal. It can be described by a power-law model, F v oc v~ T+l , with a photon index 
r = 1 .43 ± 0. 15, close to the slope Y = 1.56 ± 0.24 of the hard X-ray (E > 2.5 keV) magnetospheric component. 
The extrapolation of the X-ray magnetospheric spectrum into the optical is marginally consistent with (perhaps 
lies slightly above) the observed NIR-optical-NUV spectrum. The NUV pulsations, however, do not show a 
clear correlation with the hard X-ray pulsations. 

Subject headings: pulsars: individual (Geminga) — stars: neutron — UV: stars 
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1. INTRODUCTION 

Spin-powered pulsars show highly pulsed emission from 
the radio to 7-rays, arising from acceleration zones in their 
magnetospheres. In the UV to soft X-ray band, however, 
thermal emission from the neutron star (NS) surface can con- 
tribute significantly for middle-aged pulsars, with characteris- 
tic ages r ~ 10 4 — 10 6 yr. Spectral and timing measurements 
can separate these two components, allowing a measure of the 
surface temperature and thermal luminosity. By measuring 
thermal emission as a function of age, one can probe the equa- 
tion of state of matter at supranuclear densities in the NS core 
and constrain the surface composition. Observations with the 
Chandra and XMM-Newton X-ray observatories have begun 
to reveal much about the thermal component (see Pavlov, Za- 
vlin, & Sanwal 2002 and Kaspi, Roberts, & Harding 2004 for 
recent reviews). However, since typical effective temperatures 
of middle-aged pulsars are as low as ~ 30-100 eV, and in- 
terstellar absorption severely attenuates the flux below ~ 0.1 
keV, the X-ray observations of these objects can only probe 
the Wein tail of the surface thermal spectrum. Two issues 
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then complicate the interpretation. First, surface composition 
can dramatically affect the X-ray flux (Romani 1987; Zavlin 
& Pavlov 2002) with light element atmosphere leading to a 
large Wien excess. Second, any surface temperature inhomo- 
geneities will also complicate the spectrum, with hot spots 
disproportionally important in the high energy (X-ray) tail. 
For these reasons, comparison of the X-ray results with UV 
emission from the Rayleigh- Jeans side of the thermal bump is 
particularly valuable. The challenge here is that non-thermal 
magnetospheric emission becomes increasingly dominant as 
one moves to the red (Pavlov et al. 2002). Fortunately, the 
NUV-MAMA and FUV-MAMA detectors of the Space Tele- 
scope Imaging Spectrometer (STIS) aboard the Hubble Space 
Telescope (HST) offer access to the UV emission and provide 
the phase-resolved measurements that can help to separate the 
thermal and non-thermal fluxes. 

We report here on HST STIS observations of the middle- 
aged 7-ray pulsar Geminga. Discovered in 1972 by the SAS-2 
satellite (Fichtel et al. 1975), this object had been known only 
as a 7-ray source until it was detected in X-rays by the Ein- 
stein observatory (Bignami, Caraveo, & Lamb 1983) and as- 
sociated with a faint (V w 25.5) optical counterpart (Bignami 
et al. 1987; Bignami, Caraveo, & Paul 1988; Halpern & Tytler 
1988). The discovery of a period P = 237 ms in X-rays with 



2 



ROSAT (Halpern & Holt 1992) and 7-rays with Compton 
Gamma-ray Observatory (CSGO) (Bertch et al. 1992) proved 
the source to be a spin-powered pulsar, with a characteris- 
tic age t = P/(2P) = 342 kyr and spin-down energy loss rate 
£ = 3.3 x 10 34 ergs -1 . Contrary to most spin-powered pulsars, 
Geminga is not a strong radio source. Detection of pulsed ra- 
dio emission at 102 MHz was claimed by Malofeev & Malov 
(1997), Kuz'min & Losovskii (1997), and Shitov & Pugachev 

(1997) , but the pulsar has not been detected at other frequen- 
cies (e.g., McLaughlin et al. 1999). 

ROSAT, EUVE, and ASCA observations have established 
that the X-ray spectrum of Geminga consists of a soft ther- 
mal component, likely emitted from the NS surface, and a 
nonthermal component, presumably generated in the pulsar 
magnetosphere (Halpern & Ruderman 1993; Halpern & Wang 
1997; Jackson et al. 2002). Recent observation of Geminga 
with XMM-Newton have shown an extended emission re- 
sembling a bow-shock nebula (Caraveo et al. 2003). From 
a two-component, blackbody (BB) plus power-law (PL), fit 
of the phase-integrated XMM-Newton spectrum, Zavlin & 
Pavlov (2004a) found a temperature Tbb ~ 0.5 MK for the 
thermal component and a photon index r m 2 for the magne- 
tospheric component. The X-ray pulse profile shows a strong 
dependence on energy, changing from a single broad peak at 
E < 0.8 keV to a double-peak structure at E > 2 keV. 

The shape of the Geminga's optical spectrum remains con- 
troversial. Based on photometry with a few broad-band fil- 
ters, Bignami et al. (1996) proposed a broad emission feature 
around <~ 5000 A, superimposed on a Rayleigh-Jeans ther- 
mal spectrum (see also Mignani, Caraveo, & Bignami 1998), 
and interpreted the feature as an ion cyclotron line emitted 
from the NS atmosphere. Martin, Halpern, & Schiminovich 

(1998) reported a possible broad dip over 6300-6500 A in a 
flat (T « 1.8) spectrum spanning 3700-8000 A, but the spec- 
trum was severely contaminated by the sky background. Har- 
low, Pavlov, & Halpern (1998) detected Geminga in two near- 
IR bands, which proved that the spectrum grows toward lower 
frequencies, similar to another middle-aged pulsar B0656+14 
(Koptsevich et al. 2001). Overall, it is clear that the optical 
spectrum is predominantly nonthermal, perhaps with a hint of 
a Rayleigh-Jeans component at A < 3000 A. Optical pulsa- 
tions of Geminga were (marginally) detected in the B band 
only (Shearer et al. 1998). 

Based on three HST WFPC2 observations, Caraveo et al. 
(1996) found Geminga's parallax of 6.4 ± 1.7 mas, corre- 
sponding to d w 160 pc. Our reanalysis of these data together 
with fourth WFPC2 observation shows that the result is not re- 
liable because the exposures were too short to determine the 
Geminga's positions with the required accuracy (Pavlov et al. 
2005, in preparation). Therefore, the distance to Geminga is 
currently unknown. In this paper we will scale the distance to 
c/ = 200pc. 

Particularly interesting would be an observation of 
Geminga in the far-ultraviolet (FUV) range where one could 
expect thermal radiation from the NS surface to take over 
the apparently nonthermal radiation that prevails in the opti- 
cal. Moreover, observing pulsations of Geminga shortward of 
~ 4000 A allow one to study the transformation of the pulse 
profile in the transition from the nonthermal to thermal regime 
and elucidate the nature of the Geminga's radiation in the 
optical-UV range. To measure the spectrum and pulsations 
in the ultraviolet, we carried out an imaging observation with 
NUV-MAMA and a spectral observation with FUV-MAMA, 
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FIG. 1. — NUV-MAMA image of the field around the Geminga pulsar 
(at the center of the image). The only other point source in the field is Star 
G (e.g., Halpern, Grindlay, & Tytler 1985), used for acquisition. The inset 
shows brightness contours in the 2"1 X 2"2 region centered on Geminga. 



both with time resolution of 125 jis. These observations and 
the data analysis, including the NIR-optical-UV data, are de- 
scribed in §2. In §3 we present the spectral and timing analy- 
ses of a recent XMM-Newton observation of Geminga. Impli- 
cations of the broad-band (NIR through X-rays) observations 
are discussed in §4. The results of our work are summarized 
in §5. 

2. OBSERVATIONS AND DATA ANALYSIS 
2.1. NUV-MAMA photometry 

Geminga was observed with the STIS NUV-MAMA on 
2002 February 27 (start date is 52,332.4340 MJD UT). 
The broad-band filter F25SRF2 (pivot wavelength 2299 A, 
FWHM 1 128 A) was used in this imaging observation to min- 
imize the contribution of geocoronal lines. The data were 
taken during four consecutive orbits. The total scientific ex- 
posure time was 1 1,367 s. 

To avoid possible additional errors associated with the 
pipeline subtraction of the strong dark current background 
(see § 7.4.2 of the STIS Instrument Handbook [IHB] 1 ), we re- 
processed the "raw" NUV-MAMA images repeating all stan- 
dard calibration pipeline steps except for this subtraction. As 
an output, we obtained four flat-fielded "low-resolution" im- 
ages (1024 x 1024 pixels; plate scale 0"0244 pixel" 1 ). The 
target was detected in each of the four exposures. To increase 
the signal-to-noise ratio, S/N, we combined the images from 
the four exposures into a single image using the STSDAS 2 
task mscombine. From the sharpness of the source counts dis- 
tributions for the two point sources detected, Geminga and 
star G (see Fig. 1) we conclude that the images are aligned 
well enough for the photometry purposes. (Slight apparent 
elongations of the images of the two point sources, in differ- 
ent directions for Geminga and star G, are likely caused by 

1 STIS IHB available at http://www.stsci.edu/hst/stis/documents/handbooks/ 
currentIHB/stis_ihbTOC.html 

2 Space Telescope Data Analysis System available at 
http://www.stsci.edu/resources/software_hardware/stsdas 
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Table 1 

Numbers of NUV counts for different extraction 

APERTURE SIZES 



N, 



SN h 



N s SN S S/N e 



Fx 



(Fx) a 



1.95 1464.7 215.4 40.2 1249.2 54.5 22.9 0.496 0.222 1.26 1.33 

2.99 2119.6 502.7 21.8 1616.9 46.1 35.1 0.566 0.251 1.45 1.50 

3.95 2688.3 879.7 26.1 1808.6 50.4 35.9 0.628 0.253 1.47 1.51 

4.98 3357.4 1400.3 35.2 1957.1 57.3 34.2 0.691 0.249 1.46 1.49 

5.97 4115.1 2010.6 35.4 2104.5 58.7 35.9 0.716 0.259 1.51 1.55 

7.02 4950.4 2782.6 51.2 2167.8 70.5 30.8 0.741 0.257 1.50 1.53 

"Radius of the extraction aperture in pixels. 

b Number of background counts within the extraction aperture. 

c Source count rate corrected for the finite aperture in counts s" 1 . 

d Mean spectral flux (see eqs. [2] and [3]) in units of 10~ 18 ergs cnT 2 s" 1 



A-'. 



nonuniformities of the background.) We measure the back- 
ground, which is dominated by the detector dark current, in 
the annulus with the inner radius of 40 pixels and outer ra- 
dius of 55 pixels, centered on the source (X = 539.5 pixels, 
Y = 540.5 pixels). The mean background count rate within 
the annulus is 1.59 x 10~ 3 counts s" 1 pixel" 1 . 

To find an optimal aperture radius, we measured the num- 
ber of source counts, N s =N,-Nb (where N t is the total num- 
ber of counts, and A/j is the number of background counts 
estimated by scaling the mean background in the annulus to 
the aperture area), and its uncertainty, SN S , in apertures with 
radii of about 2, 3, 4, 5, 6, and 7 pixels (see Table 1). To 
evaluate the background uncertainty needed for calculating 
5N S , we put each of the apertures at 15 positions randomly 
distributed over the annulus, measured the number of back- 
ground counts within the aperture for each position, and cal- 
culated the root-mean-square, SNb, of the differences between 
this number and the mean background scaled to the aperture 
area. The uncertainty of the source counts can then be calcu- 
lated as SN S = [N s + 8Nl( 1 + 1 / 15)] 1/2 for each of the apertures. 
From Table 1 we see that the dependence of S/N (=N S /SN S ) 
on aperture radius has a flat maximum at a level of S/N « 35 
at r « 3-6 pixels. 

We also measured the numbers of counts in the image com- 
bined from the automatically processed images (with the dark 
current subtracted), performing standard aperture photometry 
with the IRAF task phot from the apphot package 3 . A good 
agreement with the results obtained from the direct measure- 
ments of the total (dark current plus sky) background (e.g., 
N s ± 5N S = 1821 ±51 vs. 1808 ± 50, for the 4-pixel-radius 
aperture) shows that the pipeline subtraction of the dark cur- 
rent does not introduce substantial errors in this case. 

The source spectral flux F\ is connected with the number 
of source counts in a given aperture by the integral relation 



N s = t R x XF x e x dX 



(1) 



where t is the exposure time, R\ is the integrated system 
throughput, including the Optical Telescope Assembly (OTA) 
and filter throughputs 4 , and e\ is the wavelength-dependent 

3 http://stsdas.stsci.edu/cgi-bin/gethelp.cgi?phot.hlp 

4 We corrected the throughputs supplied with the data for the time- 
dependent sensitivity loss (see http://www.stsci.edu/hst/stis/calibration/ 
reference_files/tds.html). 
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FIG. 2. — Spectral flux Fx at A = 2299 A as a function of spectral slope a\ 
in the NUV-MAMA/F25SRF2 band, for different values of E(B-V). 



encircled energy fraction. One can define the average flux in 
the filter passband as either 

Fx = ~tA „ (2) 



t J R\\e x d\ 



or 



<*a> 



te jR\\d\ ' 



(3) 



where e is in average encircled energy fraction in the filter 
passband, and N s /(te) = C is the source count rate corrected 
for the finite aperture. We calculated the average spectral 
fluxes in both ways (see Table 1) using the e\ values measured 
by Proffitt et al. (2003) for several aperture radii. We see that 
the mean fluxes, F\ ~ (F\) ~ 1.5 x 10~ 18 ergs cm" 2 s" 1 A" 1 
are close to each other for r > 3 pixels. The uncertainty of 
these values, ~ 10%, is mostly due to systematic errors in the 
encircled energy fraction. 

Another way to evaluate the flux is to assume some shape 
for the spectral flux F\ and determine its normalization mak- 
ing use of equation (1). We approximate the spectral flux 
in the F25SRF2 passband as an absorbed power law: F\ = 
F 229 9 (A/2299 ) QA io-0-4M\)E(B-v)^ where Fi229 is the intrinsic 

source spectral flux at the pivot wavelength [it coincides with 
F\ in the special case ax = 0, E(B-V)=0], and A(A) is the 
ultraviolet extinction curve (Seaton 1979). The color index 
E(B-V) is poorly known. An estimate based on the hydro- 
gen column density found from the X-ray fits (see §3.1) gives 
E(B-V)~ 0.03; below we will adopt E(B-V)=0.01-0.07 as 
a plausible range. We calculated the dependencies of F2299 
on the spectral slope ax in a reasonable range -4 < ax < 
for several values of E(B-V), based on the Af, values mea- 
sured in the 4-pixel radius aperture (see Fig. 2). We see that, 
at a given E(B-V), F2299 varies with ax by up to 20%. We 
estimate the uncertainty of the F2299 values at given ax and 
E(B-V) as w 8%-10%, mostly associated with changes of 
the MAMA imaging point spread function (PSF) between in- 
dividual observations that cause systematic uncertainties of ex 
(see Proffitt et al. 2003 and § 16. 1 of IHB). 

2.2. FUV-MAMA spectrum 

Geminga was observed with the STIS FUV-MAMA on 
2002 February 26 (start date is 52,331.2391 MJD UT). The 
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FIG. 3. — Raw FUV-MAMA spectrum of Geminga. The boxes show ap- 
proximate regions for the source and background extraction used in spectral 
analysis. 



low-resolution grating G140L (which covers the wavelength 
interval w 1150-1700 A) with the 52" x 0"5 slit was used. 
The data were taken during four consecutive orbits (including 
the target acquisition). We used a nearby field star G (V=21.3; 
Fig. 1) as the acquisition target and applied a 4 ."9 offset, de- 
duced from the positions of Geminga and the acquisition star 
measured in the archival HST images, which placed Geminga 
about 0"1 off the slit center. The total scientific exposure time 
was 10,674 s. 

For each exposure, we processed the raw "high-resolution" 
images (2048 x 2048 pixels; plate scale of 0."0122 per 
pixel — see §11 of the STIS IHB) using the calibration 
files available on 2003 July 1. As an output, we obtained 
fiat-fielded low-resolution (1024 x 1024 pixels; plate scale 
0"0244 pixel -1 ; spectral resolution 0.58 A pixel -1 ) images 
and used them for the spectral analysis. 

The processed images show a nonuniform detector back- 
ground that consists of a flat (constant) component and the so- 
called "thermal glow" component (Landsman 1998) that dom- 
inates over most of the detector area and grows with increas- 
ing the temperature of the FUV-MAMA low-voltage power 
supply (LVPS) (the average LVPS temperatures were 38.45, 
39.67, 40.89, and 41.62 C in the four consecutive orbits of 
our observation). The thermal glow is the strongest in the 
upper-left quadrant of the detector, where the dark count rate 
can exceed the nominal value, 6 x 10 -6 counts s -1 pixel -1 , by 
a factor of 20. To reduce the contamination caused by the 
thermal glow background, the source was placed close to the 
bottom edge of the detector (see Fig. 3). 

We find the Geminga's spectrum centered at Y = 105 ±2 
pixels in each of the fiat-fielded images (the centroid posi- 
tion slightly varies with X), where X and Y are the image co- 
ordinates along the dispersion and spatial axes, respectively. 
Even at this location on the detector the background still ex- 
ceeds the nominal value by a factor of 1.5-5, depending on 
the position along the dispersion axis. To improve S/N, we 
combined the images from four exposures into a single image 
using the STSDAS task mscombine. The F -positions of the 
centroids differ by less than 3 pixels for different exposures 
and different wavelengths (X -positions). 

Accurate subtraction of the enhanced, nonuniform back- 
ground [typical values are (1-3) xlO -5 counts s -1 pixel -1 ] is 
crucial to measuring the spectrum of our faint target. The 
spectral extraction algorithm implemented in the standard 
STIS pipeline (task X1D) does not adequately correct for 
the nonuniform background while extracting the spectrum of 
such a faint source and does not allow to vary the extraction 
box size with the position along the dispersion axis. There- 
fore, we used an IDL routine with additional capabilities of 
grouping and fitting the background and selecting an optimal 
extraction box size depending on the position along the dis- 
persion axis (see Kargaltsev, Pavlov, & Romani 2004). 

Since the source spectrum occupies only a small region on 



Table 2 

fuv-mama counts and fluxes in a-bins 



A- Din ( A) 


A a 


AT 

Nt 


AT 






\AI 


S/N (F x )±S(F x ) b 


1 155— 1 1&7 


9 


238.3 


157.3 


17.2 


81.1 


19.9 


4.1 


13.8 ± 3.4 


1248— 125V 


n 

y 


1 JO. / 


Oj.O 


c in 


1 

D/.l 


9.9 


6.0 


n < _i_ 1 a 
y.j ±1.4 


1 9£D 1 970 


Q 
V 


110.9 


50.7 


$ O'X 
J.jLj 


60.2 


9.5 


6.4 


7.6 ±1.2 


1271-1283 


9 


132.9 


54.0 


4.97 


78.9 


10.3 


7.7 


9.8± 1.3 


1316-1332 


7 


130.3 


57.2 


8.14 


73.1 


12.0 


6.1 


7.1 ±1.2 


1333-1347 


11 


145.8 


67.5 


8.66 


78.3 


12.6 


6.2 


8.2± 1.3 


1365-1378 


11 


96.6 


49.3 


7.11 


47.3 


10.1 


4.7 


5.3 ±1.1 


1385-1402 


11 


123.6 


62.0 


4.82 


61.6 


9.3 


6.6 


6.8 ±1.0 


1403-1431 


5 


112.5 


42.6 


6.35 


69.9 


10.6 


6.6 


6.3 ±1.0 


1432-1471 


5 


155.8 


53.1 


6.98 


102.8 


12.4 


8.3 


7.7 ±0.9 


1472-1525 


5 


140.1 


61.3 


8.01 


78.9 


12.1 


6.5 


6.0 ±0.9 


1526-1702 


5 


258.5 


155.1 


9.39 


103.4 


14.0 


7.4 


5.1 ±0.7 


Summed c 




1776.0 873.6 


28.8 


902.5 42.3 d 


20.9 


6.76±0.43 c 



"Height of extraction box, in pixels. 

b Average spectral flux and its statistical error, in units of 10~ 18 erg s -1 cm" 
A -1 , corrected for the finite aperture. 
c Values for summed A-bins. 
"Defined as E,^,,) 2 ] 1 / 2 . 

defined as [£,<F A >,AA, ± (E,(^a>?AA 2 ) 1/2 ] AA;) -1 . 



the detector, we do not attempt to subtract the background 
globally. Instead, we scan the count distribution within two 
strips, 36 < Y < 95 and 1 16 < Y < 175, adjacent to the source 
region, 96 < Y < 115. To obtain the spectrum with a suf- 
ficiently high S/N, we have to bin the spectrum heavily; af- 
ter some experimenting, we chose 12 spectral bins (A-bins; 
see Table 2). The bins exclude the regions contaminated 
by the geocoronal emission (Lya line and OI lines at 1304 
A and 1356 A) and by an artificial background structure at 
A « 1379-1384 A, Y w 96-103. The bins outside the contam- 
inated regions were chosen to have comparable S/N (w 6-8), 
whenever possible. 

For each of the A-bins, we calculate the total number of 
counts, N,, within the extraction boxes of different heights 
(one-dimensional apertures): A s = 3, 5, 7, 9, 11, 13, 15, and 
17 pixels, centered at Y = 106 for the first two A-bins and 
at Y = 105 for the rest of the A-bins. To evaluate the back- 
ground, we first clean the background strips (see above) from 
outstanding (> 10 -3 cts s -1 pixel -1 ) values ("bad pixels") by 
setting them to local average values. Then, for each of the A- 
bins, we fit the F-distribution of the background counts with 
a first-order polynomial (interpolating across the source re- 
gion), estimate the number Nt, of background counts within 
the source extraction aperture A s , and evaluate the number of 
source counts, N s = N t -Nb (Table 2). 

The uncertainty SN S of the source counts can be evaluated 
as 6N S = [N s + 5Nl(\+A s /A h )] x l 2 , where 5N b is the back- 
ground uncertainty in the source aperture. We binned the dis- 
tribution of background counts along the F-axis with the bin 
sizes equal to A s and calculated SNb as the root-mean-square 
of the differences between the actual numbers of background 
counts in the bins and those obtained from the fit to the back- 
ground. We calculated 5N S and S/N for various extraction box 
heights and found the A s values maximizing S/N for each A- 
bin (see Table 2). 

We calculated the average spectral fluxes in the A-bins (cf. 
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FIG. 4. — The measured (absorbed) FUV-MAMA spectrum of 
Geminga. The solid curve shows the best-fit absorbed blackbody model for 
E(B-V)=0.03 (T = 0.31 MK at R = 13rf 2 oo km). 
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FIG. 5. — Confidence contours (67% and 90%) for the absorbed power-law 
model fit to the FUV-MAMA spectrum, for E(B-V)=0.01, 0.03, 0.05, and 
0.07. 



eq. [3]): 



J AX R x XF x d\ _ 
J AX RxXd\ f AX Rx\d\' 



Q 



(4) 



where C, is the source count rate in the z'-th A-bin corrected 
for the finite size of the source extraction aperture, and R x is 
the system response that includes the OTA throughput and ac- 
counts for the grating and slit losses and time-dependent sen- 
sitivity losses (Bohlin, 1999; see also §3.4.12 of the HST Data 
Handbook for STIS 5 for details). The resulting flux values 
are given in Table 2, while the spectrum is shown in Figure 
4. The total flux in the 1155-1702 A range (AA = 547 A), 
can be estimated as F ~ AA (E,( /; a),AA ( ) (E, AA i) _1 - 
(3.72 ±0.24) x 10" 15 erg s" 1 cm" 2 , corresponding to the lu- 
minosity L FUV = 4ird 2 F = (1.78 ± 0.1 1) x 10 28 c/| 00 erg s" 1 . 

We fit the spectrum with the absorbed power-law model, 
Fa = ^i50o(A/1500) q * x \()-0-4M\)E(B-v)_ For plausible values 

5 http ://w w w. stsci . edu/hst/s tis/documents/handbooks/currentDHB/ 
STIS_longdhbTOC.html 
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FIG. 6. — Confidence contours (67% and 99%) for the absorbed blackbody 
model fit to the FUV-MAMA spectrum, for E(B-V)=0.03 and 0.07. 



E(B-V) = 0.01, 0.03, 0.05, and 0.07, we found the power- 
law indices a x =-3.29 ±0.53, -3.43 ±0.53, -3.56 ±0.54, 
and -3.69 ±0.54, and the normalizations F1500 = 5.38 ±0.33, 
6.25 ± 0.38, 7.26 ± 0.37 and 8.44 ± 0.52 x 10" 18 erg cm" 2 s" 1 
A" 1 , respectively (Fig. 5); the corresponding xt values are 
0.80, 0.81, 0.81, 0.82, and 0.83, for 10 degrees of freedom 
(dof). 

The inferred slope a x is close to that of the Rayleigh-Jeans 
spectrum, F x oc A" 4 , suggesting that the observed radiation is 
dominated by thermal emission from the NS surface. To es- 
timate the NS surface temperature, we fit the absorbed black- 
body model to the observed spectrum. Since the FUV fluxes 
are in the Rayleigh-Jeans part of the spectrum, the temper- 
ature is strongly correlated with the radius-to-distance ratio 
(approximately, T oc d 2 /R 2 ), as demonstrated by the confi- 
dence contours in the T-R plane (Fig. 6). For a typical NS 
radius R = 13 km and the assumed distance d = 200 pc, the 
inferred surface temperatures are 0.27 ±0.01, 0.31 ±0.01, 
0.36 ±0.02, and 0.41 ±0.02 MK, for E(B-V) = 0.01, 0.03, 
0.05, and 0.07, respectively; the corresponding x 2 values are 
0.90, 0.87, 0.85, and 0.85, for 10 dof. An example of best-fit 
blackbody spectrum is shown in Figure 4, for E(B-V)=0.03. 

2.3. NIR through FUV spectrum 

To compare the UV emission of Geminga with its NIR- 
optical emission, we plotted in Figure 7 the FUV-MAMA and 
NUV-MAMA spectral fluxes (F„) together with the fluxes at 
lower frequencies measured in eight broad passbands. Seven 
of these fluxes have been published previously (see caption to 
Fig. 7 for references), while the flux marked '555W in Fig- 
ure 7 was measured in this work from a recent observation of 
Geminga with the HST Advanced Camera for Surveys (ACS). 

Geminga was observed with the Wide Field Channel 
(WFC) of the ACS on 2003 October 7 for 6296 s total ex- 
posure (3 HST orbits, two dither positions per orbit, two ex- 
posures per dither position) in the F555W filter (ACS 'V fil- 
ter'; pivot wavelength 5358 A, FWHM 1235 A). We com- 
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FIG. 7. — NIR through FUV spectrum of Geminga. The broadband 
fluxes were measured with the HST NICMOS (F110W and F160W; Kopt- 
sevich et al. 2001), Subaru SuprimeCam (I c and R c ; Komarova et al. 2003), 
HST ACS/WFC (F555W; this work), and HST FOC (F430W, F342W, and 
F195W; Bignami et al. 1996 and Mignani et al. 1998). The solid and dash- 
dotted lines show the fits with the absorbed blackbody + power-law model 
for E(B-V)=0.03 and 0.07, respectively. The model components are shown 
by the dashed and dotted lines (see text for details). 



bined the aligned, pipeline calibrated images from the three 
orbits into a single image and performed aperture photometry 
using the phot task from the IRAF apphot package. To extract 
the source counts, we used a circular aperture with a radius of 
0."15 (3 WFC pixels) which provides an optimal S/Nw 36. 
For this aperture, the encircled energy fraction of 0.82 ± 0.04 
was determined from the empirical PSF measured for six field 
stars with apertures varying from 2 to 25 WFC pixels. The 
source count rate was corrected for the finite aperture size 
and converted to average spectral flux (cf. Eq. [3]) using the 
conversion factor (inverse sensitivity), 1 .974 x 10" 19 erg cm" 2 
A" 1 per count, for this observing mode 6 . The accuracy of the 
flux measurement, about 10%, is limited by the uncertainty 
in the encircled energy fraction and various systematic un- 
certainties. The flux we measured, (F u ) = 0.17 ±0.02 /xJy, is 
a factor of 1.6 lower than that measured by Bignami et al. 
(1996) from the HST WFPC2 observation of 1994 September 
23 with a similar filter. We remeasured the WFPC2/F555W 
flux and obtained a value consistent with our ACS result 
within the uncertainties. With the new value for the F555W 
flux and the other NIR-optical fluxes, we conclude that the 
"cyclotron feature" in the Geminga's spectrum (Bignami et 
al. 1996; Mignani et al. 1998) was likely a result of inaccurate 
photometry. 

It is obvious from Figure 7 that the NIR through FUV spec- 
trum of Geminga cannot be described by a simple power- 
law model. We fit this spectrum with a two-component, 
power-law plus blackbody, model. Since the temperature 
and the radius-to-distance ratio of the blackbody compo- 
nent are strongly correlated in the Rayleigh-Jeans regime 
(TR 2 /d 2 w const), we have to fix one of these parameters in 
the fit. For the fixed R/d = 13km/200pc, we obtained T = 
0.30 ± 0.02 MK, a v = -0.46 ± 0. 12, F = 0. 1 1 ± 0.02 /jjy for 
E(B-V)=0.03, and T = 0.41 ±0.02 MK, a„ = -0.41 ±0.13, 

6 see http://www.stsci.edu/hst/acs/documents/handbooks/DataHandbookv2 
/intro_ch34.html#l 896082 



F () = 0.12±0.02/Jy for E(B-V)=0.07 ( X l = 1.5 for 18 dof 
for each of the fits), where a v and Fq are the parameters of the 
power-law component: F v = Fq{v j\ x 10 15 Hz) a,y . Notice that 
the parameters of the blackbody component are virtually the 
same as obtained from the FUV-MAMA spectrum alone. The 
best-fit spectra and their components are shown in Figure 7. 
We see that the blackbody emission dominates at v > 1 x 10 15 
Hz (A < 3000 A), while the power-law (presumably magne- 
tospheric) emission dominates at longer wavelengths. 

2.4. Timing analysis 

For the timing analysis of the NUV-MAMA and FUV- 
MAMA data, we used the so-called TIME-TAG data files 
that contain the photon arrival times, recorded at a 125 /xs 
time resolution, and high-resolution detector coordinates (see 
§2.2) associated with each of the events. We use 2688 
NUV-MAMA events extracted from an aperture of 8 high- 
resolution pixels radius (includes 66% of source counts) and 
1939 FUV-MAMA events extracted from the above-defined 
A-bins with heights of extraction box varying from 14 to 22 
high-resolution pixels, depending on the A-bin (includes 46% 
of source counts). The arrival times are corrected for the 
Earth and spacecraft motions and transformed to barycen- 
tric dynamical times (TDB) at the solar system barycenter, 
using the STSDAS task odelaytime. The time spans of the 
FUV-MAMA and NUV-MAMA observations are 19,181 s 
and 19,785 s, respectively, with a gap of 84,041 s between 
the last FUV-MAMA event and the first NUV-MAMA event. 

The expected frequency of Geminga's pulsations at the 
epoch of our observation, around 52,332 MJD, can be esti- 
mated from the previous timing observations in 7-rays and 
X-rays. The most recent ephemerides of Geminga were 
published by Jackson et al. (2002; J02 hereafter). These 
authors found a small glitch in the Geminga's timing his- 
tory and presented a post-glitch ephemeris for a time inter- 
val of 50,382-51,673 MJD. Although our observation was 
taken 659 days after the end of that interval, extrapolation 
of this ephemeris to 52,332 MJD predicts the frequency, 
/j02 = 4,217,608.6953/iHz, with a formal uncertainty of 
±0.0013 /iHz, that is about three orders of magnitude smaller 
than we can achieve in our relatively short observation (see 
below). Therefore, we adopt /jo2 as an estimate of ex- 
pected frequency and look for pulsations in its vicinity, in a 
/j 2 ± (2r span ) _1 frequency range. 

Since the longer time span, r span = 123,005 s, of the joint 
FUV+NUV data set allows a tighter constraint on the pul- 
sation frequency, we start from the analysis of this data set. 
First, we apply the Z 2 test (Buccheri et al. 1983), calculat- 
ing the Z 2 statistic as a function of trial frequency in the 
range of /}o2 ±4 /iHz for n = 1 through 8, where n is the 
number of harmonics included. For each of the n values 
examined, we found statistically significant pulsations, with 
frequencies of Z 2 maxima within (-0.3,+0.7)^Hz around 
/j02. The most significant result is obtained for n = 6: max- 
imum Z 2 = 53.1 at / = 4,217, 608.8/iHz; the probability 
to obtain this value by chance is 4 x 10~ 7 . To better es- 
timate the uncertainty of pulsation frequency, we also ap- 
plied the odds-ratio method of Gregory & Loredo (1992; 
see also Zavlin, Pavlov, & Sanwal 2004) and found / = 
4,217,608.3 (-0.9, +0.7) /zHz for the median frequency and 
68% uncertainties, and / = 4,217,608.8 ± 1.2/iHz for the 
mean frequency and standard deviation. Within the uncertain- 
ties, these frequencies virtually coincide with the frequency 
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FIG. 8. — UV light curves of Geminga folded with the J02 ephemeris. Top: 
NUV-MAMA light curve, obtained using the data from all four orbits. The 
estimated average background level in the 20-bin light curve is 44.3 counts 
per phase bin. Middle: FUV-MAMA light curve, obtained using the data 
from all four orbits. Bottom: FUV-MAMA light curve, obtained using only 
the data from the first two orbits. The dashed lines with associated error bars 
show the corresponding background light curves (see text). 



predicted by the J02 ephemeris. The expected frequency shift 
during the FUV+NUV observation, fT sp . dn = 0.024 ^Hz, is 
much smaller than the frequency uncertainties in our mea- 
surement, which means that this observation is not sensitive 
to frequency derivative /. 

Since our spectral analysis has shown that the FUV radia- 
tion is predominantly thermal, while the NUV radiation has 
a significant contribution from the magnetospheric compo- 
nent (see §2.3 and Fig. 7), one can expect different strength 
and shape of pulsations in the FUV and NUV bands. There- 
fore, we analyzed these two data sets separately. Since the 
frequency we measured from the FUV+NUV data is consis- 
tent with the J02 ephemeris, and an XMM-Newton obser- 
vation taken 37 days later also suggests that the ephemeris 
may still be valid (see §3.2), we folded the times of arrival 
with the J02 ephemeris, choosing the same zero-phase epoch, 
50,382.999999364 MJD. 



The folded (source plus background) light curve in the 
NUV-MAMA band, plotted in the upper panel of Figure 8, 
shows one broad (FWHM « 0.8 in phase), flat-top peak per 
period, centered at <j> w 1 .0. The most notable feature of the 
pulse profile is the narrow dip at (j> « 0.45. The Z\ test shows 
that the pulsations are statistically quite significant, with the 
main contribution coming from the fundamental frequency: 
Z\ = 22.3 corresponds to 1 x 10~ 5 probability of false result. 
The pulsed fraction, defined as the ratio of the number of 
counts above the minimum level to the total number of counts 
in the light curve, is about 28%, which corresponds to the in- 
trinsic source pulsed fraction / p w 40%. 

The Zl test for the FUV-MAMA data set shows most signif- 
icant pulsations for n = 4: Zj = 25.08 corresponds to 99.84% 
(3.2 a) significance. The lower significance of the FUV pul- 
sations, compared to the NUV pulsations, can be caused by 
a lower S/N in the spectroscopic mode. Since the thermal- 
glow background was growing with increasing the LVPS tem- 
perature in the course of our observation (see §2.2), we per- 
formed the timing analysis for various combinations of orbits 
and found that indeed the pulsations were more significant in 
earlier orbits. For instance, in the first two orbits (r span = 7638 
s,N = 844 counts) the most significant Z\ = 23 .46 corresponds 
to 99.990% (3.9cr) significance. 

The four-orbit and two-orbit FUV-MAMA light curves are 
shown in the middle and lower panels of Figure 8. In the same 
panels we show the light curves for the background counts 
extracted from two boxes centered at Y = 180 and Y = 240 
high-resolution pixels, with the same heights as used for ex- 
traction of the source events. The background light curves 
do not show statistically significant pulsations. Both the four- 
orbit and two-orbit source-plus-background light curves show 
a sharp, asymmetric dip at approximately the same phase as 
the NUV-MAMA light curve. A hint of a shallower dip, better 
pronounced in the two-orbit light curve, is seen at 4> w 0.95. 
The pulsed fraction in the observed (source + background) 
radiation is about 35% and 45% for the four orbits and first 
two orbits, respectively. The corresponding intrinsic pulsed 
fractions are rather high, about 60%-70%. It should be noted, 
however, that these values are rather uncertain because of the 
large statistical error of the minimum level. 

3. X-RAY SPECTRUM AND PULSATIONS OF GEMINGA 

To better understand the UV spectrum and pulsations of 
Geminga, observations at X-ray wavelengths are particularly 
useful. The deepest observation of Geminga in X-rays was 
carried out with the XMM-Newton observatory on 2002 April 
4-5 (orbit 425). The EPIC 7 -MOS and EPIC-pn instruments 
observed the pulsar for 101.4 and 71.4 ks of effective expo- 
sures, respectively. Two EPIC-MOS detectors were operated 
with medium filters in Full Frame mode providing an image 
of a large area, r ~ 14', with time resolution of 2.6 s. EPIC- 
pn was used in combination with thin filter in Small Win- 
dow mode which covers a 4.'37 x 4f37 region and provides 
a 5.7 ms time resolution. First results of this observation have 
been reported by Zavlin & Pavlov (2004a) and Caraveo et al. 
(2004a,b). Here we briefly describe the X-ray spectrum and 
pulsations of Geminga, with emphasis on the properties most 
useful for the comparison with the optical-UV data. 

3.1. X-ray spectrum 
7 European Photon Imaging Camera. 
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FIG. 9. — Confidence contours (68%, 90%, and 99%) in the temperature- 
radius plane obtained from fitting the EPIC-pn spectra (solid lines) with the 
TS+TH+PL and BB+PL models (labels near the contours). The TS+TH+PL 
contours were obtained at the fixed parameters of the PL component; nn was 
free for the upper contours, while it was fixed at the value obtained from 
ROSAT data for the lower contours. All the model parameters were free for 
the EPIC-pn BB+PL contours. The dashed and dash-dotted lines show the 
confidence contours obtained from fitting the FUV-MAMA spectrum with a 
blackbody model for two values of the color index E(B-V). 



The most detailed X-ray spectrum of Geminga was ob- 
tained with the EPIC-pn instrument. The EPIC-pn data, pro- 
cessed with the SAS package 8 (ver. 6.0.0), were used for the 
spectral and timing analysis. We extracted the source (plus 
background) photons from a 40"-radius circle centered at the 
pulsar position, which contains about 88% of source counts. 
The estimated total source count rate (corrected for finite ex- 
traction radius) is 0.813 ±0.004 counts s" 1 in the 0.2-10 
keV range for single and double events (with photon-induced 
charge detected in a single CCD pixel and two adjacent pix- 
els). The 0.2-10 keV phase-integrated spectrum was binned 
in 222 spectral bins with at least 40 source counts per bin. 
The detector response matrix and effective area were gen- 
erated with the rmfgen and arfgen tools, respectively. Fit- 
ting this spectrum with a two-component, blackbody (BB) + 
power-law (PL), model, we find the blackbody temperature 
r bb = 0.47 ± 0.02 MK and radius R = (17.0 ± 2.5) d 200 km, 
which suggests that the thermal component originates from 
the NS surface. The PL component, with a photon index T = 
2.02 ± 0.05 dominates at energies E > 0.6 keV and contains 
about 10% of the total luminosity in the 0.2-10 keV band, 
rP 1 



^0.2-10keV 



(2.6 ±0.1) x lO^fooergss" 1 ~ 7 x lO^dj^E 
Extrapolated into the optical domain, the PL component ex- 
ceeds the observed optical fluxes by a factor of 100-500, 
which might be interpreted as a flattening of the pulsar mag- 
netospheric spectrum at lower photon energies. The hydrogen 
column density derived from this fit is n H = (2.9 ±0.2) x 10 20 
cm L . 

Although the two-component model cannot be rejected 
based on the overall fit quality (x 2 = 1.11 for 217 dof; sys- 
tematic errors in the EPIC-pn response not included), the fit 
residuals show some excess of observed counts over the best- 
fit model at higher energies, E > 7 keV, indicating a harder PL 
spectrum. Indeed, fitting the high-energy tail (E > 2.5 keV) 

8 http : / /xmm . vil spa. esa.es 



(2.2 ±0.2) x 10 30 c/| 00 erg s" 1 (xl = 0.98 for 30 
dof). The BB+PL fit with T fixed at this value is statisti- 
cally unacceptable {\ v = 2.95 for 219 dof). Therefore, we 
tried a three-component model consisting of soft (TS) and 
hard (TH) blackbody components and a PL component. With 
the PL parameters fixed at the values obtained from the best 
PL fit in the 2.5-10 keV band (V = 1.56, N = 5.5 x 10" 5 
photons cm" 2 s" 1 kev" 1 at 1 keV), we obtain the follow- 
ing parameters for the thermal components: T s = 0.49 ±0.01 
MK, R s = (12.9 ± 1.0)^200 km, T h = 2.32 ± 0.08 MK, R h = 
(46 ± 12)c/ 2 oo m, and « H = (2.4 ±0.2) x 10 20 cm" 2 (x 2 =1.10 
for 217 dof). Fitting the EPIC-MOS spectra obtained in this 
observations yields almost the same model parameters (dis- 
carding the EPIC-MOS events below 0.3 keV, where the re- 
sponses of the MOS detectors are known very poorly). 

In the TS+TH+PL model the TS component can be inter- 
preted as emission from the bulk of NS surface, the TH com- 
ponent can be ascribed to emission from smaller, hotter re- 
gions of the NS surface, and the PL component represents 
the magnetospheric radiation. Such an interpretation of the 
Geminga's X-ray spectrum is in line with the results obtained 
from Chandra and XMM-Newton observations of the other 
bright middle-aged pulsars, B0656+14 and B1055-52, whose 
X-ray spectra can also be desrcibed by the TS+TH+PL model 
with similar parameters (Pavlov et al. 2002; Zavlin & Pavlov 
2004a). However, the effective radius of the Geminga's TH 
component is much smaller than those of B0656+14 and 
B1055-52, about 0.6(d/0.3kpc)km and 0.4(c//0.7kpc)km, 
respectively. 

The confidence contours for the temperature and radius of 
the TS component are shown in Figure 9. In the same fig- 
ure we plotted the temperature-radius confidence contours ob- 
tained from the blackbody fit of the FUV-MAMA spectrum, 
for E(B-V)=0.03 and 0.07 (see §2.2 and Fig. 6). We see that 
at plausible values of interstellar extinction, E(B-V)< 0.07, 
the FUV contours lie at smaller radii (or lower temperatures) 
than the X-ray contours. This means that the extrapolation of 
the thermal X-ray component into the UV-optical goes above 
the observed FUV flux: (TR 2 ) X /(TR 2 ) FUV ~ 1.8, 1.6, 1.4, 
and 1.2 for E(B-V) = 0.01, 0.03, 0.05, and 0.07, respectively, 
for the best-fit parameter values (see also Fig. 10). Such be- 
havior is in contrast to other neutron stars observed in both 
X-rays and optical, for which such an extrapolation usually 
underpredicts the optical-UV fluxes 9 by a factor of 2-7 (see 
§4.1). To obtain a similar ratio for Geminga, one would have 
to assume unrealistically high extinction, E(B-V) ~ 0.2-0.4. 

If we adopt the above-described two-component (BB+PL) 
model, the discrepancy between the X-ray and FUV temper- 
ature and radius is even more pronounced, as demonstrated 
by the corresponding confidence contours in the upper left 
part of Figure 9 and the upper panel of Figure 10. We note, 
however, that the «h values corresponding to these contours, 
(2.4 ±0.2) and (2.9 ±0.2) x 10 20 cm" 2 for the TS+TH+PL 
and BB+PL models, respectively, significantly exceed the 
n H ^ (1.1 ±0.2) x 10 20 cm- 2 obtained from the ROSATPSPC 
observations (Halpern & Wang 1997), which indicates a dis- 
crepancy between the PSPC and EPIC responses at low ener- 
gies. If we fix n H at the best-fit PSPC value, the confidence 
contours shift to higher temperatures and lower radii (see the 
lower-right EPIC-pn contours in Fig. 9), overlapping the FUV 

9 A possible exception is the Vela pulsar (see Romani, Kargaltsev, & 
Pavlov 2005). 
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FIG. 10. — NIR through X-ray spectrum of Geminga for different X-ray 
spectral models and different color indices. The solid lines show the best-fit 
(absorbed) spectra in the X-ray range and their extrapolations into the NIR- 
FUV range. The short-dash and dash-dot lines show the (soft) thermal and 
PL components, respectively, the dotted lines in three lower panels show the 
TH component (its contribution is negligible in the NIR-FUV range), and 
the long-dash lines present the unabsorbed total spectra. The crosses depict 
the measured NIR-FUV spectral fluxes (cf. Fig. 7). The hatched areas along 
the PL and thermal components in the NIR-FUV range demonstrate prop- 
agated uncertainties of the corresponding extrapolations. The upper panel 
shows a two-component (BB+PL) X-ray fit, while three lower panels show 
TS+TH+PL fits with fixed parameters of the PL component. The fit shown 
in the lower panel was obtained at a fixed «h value, while «h was a fitting 
parameter in three upper panels. (See text for more details.) 



contours. For this tin, the FUV fluxes lie on the extrapola- 
tion of the best-fit X-ray TS component at E(B-V) < 0.04 
(see lower panel of Fig. 10). Since neither EPIC nor PSPC 
have been accurately calibrated for very soft spectra, system- 
atic errors can substantially exceed the statistical errors, and 
the model parameters inferred from such fits may not be very 
accurate. Therefore, there still remains some uncertainty in 
the comparison of the X-ray (Wien) and UV (Rayleigh-Jeans) 
tails of the thermal spectrum. However, even with account for 
this uncertainty, Geminga exhibits a fainter UV-optical ther- 
mal radiation, relative to the soft X-ray radiation, than the 
other neutron stars for which such a comparison is possible 
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FIG. 11 .— ASCA GIS and XMM EPIC-pn light curves in the 0.5-2 keV 
range folded with the J02 ephemeris. 



(see §4.1). 

Figure 10 also shows that the continuation of the best-fit 
X-ray PL into the optical very strongly overpredicts the ob- 
served NIR-optical fluxes for the BB+PL model. However, 
the predicted and observed fluxes become marginally consis- 
tent if we use the PL component inferred from the E > 2.5 
keV spectral tail. 

3.2. X-ray pulsations 

To study the X-ray pulsations of Geminga, we use the same 
EPIC-pn data (r span = 101.9 ks, epoch of the middle of the 
time span 52,369.2997 MJD). First, we measured the pulsa- 
tion frequency using the and odds-ratio methods, for vari- 
ous energy bands and extraction radii, and found most prob- 
able frequencies in the range of 4,217,607.75^,217,607.96 
/iHz, with typical uncertainties of about 0.1 /iHz for indi- 
vidual measurements. For example, the odds-ratio method 
applied for 42, 170 events in the 0.23^1.0 keV band, extracted 
from 40"-radius circle, gives / = 4,217,607.85 ± O.lO^Hz 
for the mean frequency and standard deviation, and / = 
4,217, 607.86(-0.05,+0.05;-0.16,+0.12;-0.24,+0.19)/zHz 
for the median frequency and 68%, 90%, and 99% un- 
certainties. The most probable frequencies are con- 
sistently lower, by 0.1-0.3 /xHz, than the frequency 
/ J02 = 4,217,608.0664 ±0.0013^Hz predicted by the J02 
ephemeris. However, since the differences do not exceed 
3 a uncertainties of our measurements, it is still possible 
that the J02 ephemeris is applicable at the epoch of the 
XMM-Newton observation. 

We have also directly checked the phase alignment of the 
light curves extracted from the XMM-Newton data with those 
observed by ASCA in 1999 October 5-11. The lower panel 
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of Figure 1 1 shows the XMM-Newton light curve folded with 
the J02 ephemeris, in the energy band 0.5-2 keV (10,264 
counts in 40"-radius aperture). The upper panel of the same 
figure shows the light curve obtained with the two ASCA GIS 
instruments (r exp = 207.8 ks, r span = 486.5 ks, epoch of the 
middle of the time span MJD 51,459.7356; 1,819 counts in 
3'-radius aperture) folded with the same ephemeris and in the 
same energy range. We see that not only the shapes of these 
light curves are virtually the same, but also their phases are in 
excellent agreement, within the phase uncertainty (~ 0.12) of 
the J02 timing solution propagated to the epoch of the XMM- 
Newton observation. Therefore, we will assume that the J02 
ephemeris is still applicable in 2002 April and use it to com- 
pare the light curves observed with different instruments. 

The background-subtracted light curves in the energy 
ranges 0.2-0.5, 0.6-1.0, and 2-8 keV are shown in Figure 
12. In the 0.2-0.5 keV and 2-8 kev bands the radiation is 
dominated by the TS and PL components, respectively, while 
the 0.6-1 .0 keV band was chosen around the maximum of the 
TH component (see Fig. 10). The light curves were extracted 
from a smaller 30"-radius aperture and a shorter, 80.0 ks, time 
span (exluding intervals of strong background flares at the be- 
ginning and end of the observation) to reduce the background 
contamination and maximize the signal-to-noise ratio. The 2- 
8 keV light curve (pulsed fraction / p = 34% ±8%) shows two 
pronounced peaks per period, resembling the 7-ray light curve 
(albeit with smaller distance between the peaks) and a hint of 
a third peak, at <\> w 0.2. On the contrary, the 0.2-0.5 keV light 
curve (/ p = 30% ± 2%) is characterized by one broad peak per 
period (with small "ripples", perhaps due to contribution from 
the PL and TH components). The 0.6-1.0 keV light curve 
shows the highest pulsed fraction, / p = 62% ± 5%, with one 
asymmetric peak, possibly comprised of several peaks asso- 
ciated with contributions from different components in this 
band (Fig. 10). The minimum of the 0.2-0.5 keV light curve 
is approximately aligned in phase with one of the minima of 
the 2-8 keV light curve, being shifted by A<j> w 0.1 from the 
sharp dips of the NUV and FUV light curves. (One should 
remember, however, that the shift can be caused by errors in 
phase alignment.) Examples of X-ray light curves for other 
energy ranges can be found in Zavlin & Pavlov (2004a) 10 and 
Caraveo et al. (2004a,b). A more detailed discussion of the 
thermal and nonthermal light curves is presented in §4. 

4. DISCUSSION 

The above-described results of the observations of Geminga 
show that both the NIR-optical-UV and X-ray emission are 
comprised of thermal and nonthermal components, with quite 
different spectra and light curves. In the following, we dis- 
cuss the multiwavelength properties of these components sep- 
arately. 

4.1. Thermal component(s) of the Geminga 's emission 
4.1.1. Spectrum 

It follows from §2.2 and §3.1 that the Geminga's radiation 
is predominantly thermal at 4eV < E < 0.5 keV. The ob- 
served FUV and soft X-ray (TS component) spectra represent 
the Rayleigh-Jeans and Wien tails of the thermal spectrum 
emitted from the NS surface. The blackbody fits of the soft 
X-ray emission give the NS surface temperature in a range 

10 The inaccurate estimates of pulsed fractions given in Figure 8 of Zavlin 
& Pavlov (2004a) should be disregarded. 
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FIG. 12. — Background-subtracted light curves of Geminga in UV (NUV- 
MAMA + FUV-MAMA), X-rays (EPIC-pn) and 7-rays (EGRET) bands, 
folded with the J02 ephemeris. The 7-ray light curve is taken from Jackson 
et al. (2002). 
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FIG. 13. — Temperature-radius confidence contours (68%, 90%, and 99%) 
for the isolated neutron star RX J 1856.5-3754 obtained from the X-ray ob- 
servations with different instruments (solid contours) and FUV-MAMA ob- 
servations (dash and dash-dot lines). 



of 0.45-0.53 MK (« 39^6 eV), corresponding to effective 
radii of 21-8 km, at d = 200 pc (see Figs. 9 and 10). The 
uncertainty in these parameters is mostly due to the poorly 
calibrated responses of the EPIC detectors at low energies. 
Moreover, these temperatures are somewhat lower, and the 
radii larger, than those estimated from the previous ROSAT 
and EUVE observations (Halpern, Martin, & Marshall 1996; 
Halpern & Wang 1997). Since significant variations of the 
NS temperature and emitting area in a ten year span of these 
observations can hardly be expected, the discrepancy is most 
likely due to discrepant instrument responses at low energies. 

As we have shown in §3.1, the observed thermal UV spec- 
trum of Geminga either matches the continuation of the ther- 
mal X-ray spectrum or lies somewhat below that continuation 
(up to about one stellar magnitude), depending on assumed 
extinction and X-ray spectral model. On the contrary, other 
NSs observed in both UV and optical show relatively brighter 
Rayleigh-Jeans components, well above the continuation of 
the X-ray thermal spectrum. To demonstrate this difference 
and show that it is not associated with uncertainties in instru- 
ment responses, we re-analyzed the FUV-MAMA and EPIC 
data on the best-studied isolated neutron star, RX J 1856.5- 
3754 (J1856 hereafter; see Trumper et al. 2004 for a recent 
review of its properties). We used the FUV-MAMA obser- 
vation of 2002 October 26 (exposure time 13,451 s), ana- 
lyzed the data as described in §2.2, and confirmed that the 
spectrum follows a Rayleigh-Jeans law (Pons et al. 2002), 
with a total flux F = (1.89 ±0.09) x 10" 14 ergs s" 1 cm" 2 
in the 1155-1702 A range. Fitting the FUV-MAMA spec- 
trum with a blackbody model gives T = (0.45 ± 0.02) /^if^o 
MK and (0.55 ±0.02) R\jdj 20 MK, at plausible color indices 
E(B-V) = 0.01 and 0.03 (R u = fl/13km, d m = t//120pc). 
The temperature-radius confidence contours of these fits are 
shown in Figure 13. We also re-analyzed the archival XMM- 
Newton observation of 2002 April 8-9 together with the re- 
cent observation of 2004 April 17-18. In the observations of 
2002 and 2004, the EPIC-pn was operated in Small Window 
mode with thin filter (40.0 ks effective exposure) and Tim- 
ing mode with thin filter (64. 1 ks effective exposure), respec- 



tively. The EPIC-MOS observations of 2004 were carried out 
in Full Frame mode with thin filter (the same 65.3 ks effec- 
tive exposures for MOS1 and MOS2). We did not use the 
EPIC-MOS observations of 2002 because MOS1 was oper- 
ated in Timing mode, which is very poorly calibrated for this 
instrument, and MOS2 was operated in Small Window mode, 
with a field-of-view 100" x 100" (for the central CCD) too 
small to reliably subtract the background. We found that the 
two EPIC-pn spectra of J 1856 are well consistent with each 
other in the 0.3-1 keV range (there is no spectral informa- 
tion available below 0.3 keV in the data collected in Timing 
mode). Since J 1856 does not show a nonthermal component, 
we fit the spectra with a single-component BB model and plot 
the corresponding confidence contours in Figure 1 3 [«h 20 = 
0.66±0.3,0.04(-0.04,+0.12),and0.03(-0.02,+0.41)forthe 
pn, MOS1 and MOS2 detectors, respectively]. In the same 
figure we also plot the confidence contours obtained from fit- 
ting the 449.9 ks observation of 2001 October 8-15 with the 
Low Energy Transmission Grating Spectrometer (LETGS) on 
Chandra (see Burwitz et al. 2003); the corresponding hydro- 
gen column density is n H ,20 = 0.86 ± 0.15. We see that the 
observations with different X-ray instruments yield quite dif- 
ferent spectral parameters. This demonstrates once more the 
lack of proper cross-calibration of instrument responses to 
soft spectra and the fact that systematic uncertainties greatly 
exceed statistical ones for spectra with good statistics. 

Even with allowance for the uncertainties in instrument re- 
sponses, we see from Figure 13 that the UV contours lie well 
above the X-ray contours, i.e., the extrapolation of the X- 
ray blackbody spectrum of J 1856 into the UV-optical range 
strongly underpredicts the observed UV-optical fluxes, con- 
trary to Geminga (cf. Fig. 9). The UV-optical excess in the 
thermal spectra of J 1856 and other so-called X-ray Dim Iso- 
lated NSs (e.g., RX J0720.4-3125) could be explained as- 
suming that X-rays are emitted from a small hotter area while 
the optical-UV radiation is emitted from the bulk of NS sur- 
face, including colder areas invisible in X-rays (e.g., Pavlov 
et al. 2002). Obviously, the apparently smaller UV-emitting 
area of Geminga, as compared to the X-ray-emitting area, 
cannot be explained by a nonuniform temperature distribu- 
tion. We might speculate that the temperature distribution 
over the bulk of Geminga's surface is more uniform than in 
the case of J 1856, e.g., because of a different geometry and 
strength of the magnetic field that affects the heat conductiv- 
ity and, hence, the surface temperature distribution. However, 
to explain why the more uniformly heated Geminga exhibits 
quite substantial pulsations of its thermal radiation (/ p « 30% 
in soft X-rays) while no pulsations have been detected from 
J 1856, one has to assume a special orientation of the spin axis 
ofJ1856. 

Any realistic interpretation of thermal emission from NSs 
should take into account possible deviations of thermal spec- 
tra emitted from NS surface layers (e.g., atmospheres) from 
the idealized Planck spectra as well as the anisotropy of the 
surface emission associated with strong magnetic fields (e.g., 
Pavlov et al. 1995; Rajagopal, Romani, & Miller 1997; Zavlin 
& Pavlov 2002). For instance, since the X-ray spectrum emit- 
ted from a strongly ionized hydrogen atmosphere is harder 
than the Planck spectrum, a blackbody fit of such a spectrum 
gives a temperature exceeding the actual effective temperature 
by a factor of 1.5-2.5 and a radius a factor of 3-15 smaller 
than the actual radius of the NS, at a given distance. More- 
over, the optical part of the spectrum emitted from such an 
atmosphere strongly exceeds the extrapolation of the black- 
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body fit of its X-ray spectrum into the optical domain (Pavlov 
et al. 1996). Although the fully ionized atmosphere models 
are not applicable to cold NSs (e.g., they give an improbably 
large radius-to-distance ratio, R ~ 100 km at d = 200 pc, for 
Geminga; Meyer, Pavlov, & Meszaros 1994), atmospheric ef- 
fects might, in principle, explain the large difference between 
the observed UV spectrum and the continuation of the black- 
body fit of the X-ray thermal spectrum in J 1856 and similar 
NSs. However, no realistic models adequately describing the 
observed broad-band spectrum of J 1856 have been suggested 
so far, which is not suprising given the extremely complicated 
physics of the dense, strongly magnetized matter at the rel- 
atively low temperature of the surface layers. If we adopt 
such an interpretation of the strong deviation of the J 1856 
broad-band spectrum from a pure blackbody spectrum, then 
we have to explain why the broad-band spectrum of Geminga 
is so different from that of J 1856. Possible hypotheses might 
involve different chemical compositions of the surface lay- 
ers and/or substantially different magnetic fields (a crude es- 
timate of the Geminga's magnetic field is ~ 2 x 10 12 G, but 
the magnetic field of J1856 is quite uncertain 11 ). Moreover, 
the surface layers of Geminga and J 1856 might be in differ- 
ent phase states. For instance, one could speculate that the 
cold surface of Geminga is in a solid state while the (hotter) 
surface of J 1856 is in a gaseous or liquid state, which might 
explain their different spectra. To distinguish between these 
possibilities, reliable models for NS thermal emission at rela- 
tively low temperatures, with account for the contribution of 
molecules in the opacity of gaseous atmospheres with strong 
magnetic fields (Turbiner & Lopez Vieyra 2004) and possi- 
ble condensation of the surface layers into a liquid or solid 
state (van Adelsberg, Lai, & Potekhin 2004), are to be devel- 
oped and compared with the observational data. Until reli- 
able models are available, the temperatures and radii obtained 
from applying simplified models (blackbody, fully ionized at- 
mospheres, partially ionized atmospheres without molecules) 
should be considered as crude estimates only, and any con- 
clusions based on such fits should be considered with caution. 
However, although we cannot trust absolute values of the pa- 
rameters obtained with the aid of simplified models, some in- 
teresting qualitative results can be obtained from comparison 
of the same parameter measured for different NSs. For in- 
stance, fits of the soft X-ray spectra of Geminga and an older 
pulsar B 1055-52 with any model available give a lower tem- 
perature for the younger Geminga, which may have very in- 
teresting implications for the NS cooling models, suggesting 
different masses of these NSs (Yakovlev & Pethick 2004). 

As we mentioned in §3.1, in addition to the thermal soft 
(TS) component, the X-ray spectrum of Geminga apparently 
has a thermal hard (TH) component, with a much higher tem- 
perature, Th~2 MK, and an apparent (isotropic) luminosity 
La~4 x 10 29 ergs s" 1 . Although such a component has been 
seen in the spectra of other middle-aged pulsars, the effective 
radius, Rh ~ 50c/200 m, for the Geminga's TH component is 
surprisingly small in comparison with the conventional po- 
lar cap radius, fl pc = (27r/7? 3 ) 1/2 c~ 1/2 w 300 m, suggested by 
the pulsar models. Such a small value of Rh might be ex- 
plained by a projection effect (if the magnetic axis remains 
almost perpendicular to the line of sight in the course of NS 
rotation), but this explanation can hardly be reconciled with 

11 For instance, Pavlov & Zavlin (2003) consider a possibility that J 1856 
is a millisecond pulsar with a very low magnetic field, B ~ 10 s — 10 9 G, while 
Trumper et al. (2004) suggest that it has a very strong field B > 10 13 G. 



the high pulsed fraction at energies where the TH component 
contribution is maximal (see §4.1.2). On the other hand, we 
should remember that the TS component was obtained assum- 
ing Planck spectra for both thermal components and a single 
power law for the magnetospheric spectrum. Because both 
these assumptions are not necessarily correct, we cannot rule 
out the possibility that the "TH component" is simply asso- 
ciated with a harder high-energy tail of the surface radiation 
(compared to the pure Wien spectrum) or it is due to a steepen- 
ing of the slope of the phase-integrated magnetospheric spec- 
trum with decreasing photon energy (see §4.2). 

4.1.2. Pulsations in thermal emission 

One of the most intriguing results of our STIS MAMA ob- 
servations of Geminga is the strong, non-sinusoidal pulsations 
in the FUV range, where the spectrum is dominated by the 
thermal component, most likely emitted from the bulk of NS 
surface. The shape of the FUV pulsations is different from 
that of the soft X-ray pulsations, where the TS component 
dominates (see Fig. 12). Obviously, neither FUV nor soft X- 
ray pulsations can be produced by the locally isotropic black- 
body emission. To explain the unusual pulse shape and the 
large pulsed fraction of the thermal FUV and soft X-ray ra- 
diation, we have to invoke effects of strong magnetic field 
on the angular dependence of NS surface emission or assume 
that there is a "screen" in the NS magnetosphere which may 
partially eclipse the surface emission at some rotation phases. 

In a strong magnetic field, B> \0 n (E/\ keV) G, when the 
electron cyclotron energy E c exceeds the photon energy, the 
local emission is essentially anisotropic (in particular, beamed 
along the direction of the magnetic field), which may lead to 
strong pulsations of the thermal radiation. The angular dis- 
tribution and the shape of pulsations depend on the proper- 
ties of the emitting region. For instance, the angular distri- 
bution of local emission from a fully ionized NS atmosphere 
shows a strong, narrow peak [A6 ~ (E/E c )^ 2 ] along the mag- 
netic field (pencil component) and a broad fan-like component 
across the magnetic field (Pavlov et al. 1994). When inte- 
grated over the visible surface of a NS with a dipole magnetic 
field, the angular distribution of NS radiation is beamed along 
the magnetic axis, even in the case of a uniformly heated NS 
surface (Zavlin & Pavlov 2002). Such peaks could explain the 
soft X-ray pulsations (at E > kTes), including the observed in- 
crease of pulsed fraction with energy. In this hypothesis, the 
0.2-0.5 keV pulse profile (see Fig. 12) can be interpreted as a 
sum of a smooth thermal component (the broad thermal peak, 
with a maximum at <f> « 0.8, corresponds to the closest ap- 
proach of the magnetic axis to the line of sight), and small 
"wiggles" (e.g., at <f> w 0.9) due to the contribution of the 
magnetospheric radiation. However, at E <C kT e g the peaks 
in the model angular distribution are too low to be respon- 
sible for the observed FUV pulsations. On the other hand, 
as we mentioned above, the fully ionized atmosphere mod- 
els are not directly applicable to the cold Geminga, while the 
partially ionized atmospheres have not been well investigated. 

If the NS surface matter is in a condensed state, we also 
should expect an anisotropic emission. Although the angu- 
lar distribution of emission from a condensed surface has not 
been studied, the examples of spectral emissivity for several 
directions, calculated by van Adelsberg et al. (2004), sug- 
gest that at least local radiation is beamed along the magnetic 
field. To understand whether the radiation from the entire NS 
surface can show pulsations similar to those observed from 
Geminga, the local specific fluxes should be integrated over 
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the visible NS surface, for various magnetic field geometries 
and orientations of the spin and magnetic axes (see Perez- 
Azorin, Miralles, & Pons 2004 for a few examples). 

An alternative explanation for the narrow deep minima in 
the UV pulse profiles could be a partial eclipse by an ob- 
ject co-rotating with the NS. Since the shapes of the UV and 
soft X-ray light curves are different (in particular, the min- 
ima are broader in soft X-rays), the eclipsing object should 
have a wavelength-dependent effective size. A natural candi- 
date for such a screen is the magnetospheric electron-positron 
plasma which can absorb the NS surface radiation as a re- 
sult of the cyclotron resonance scattering in a resonance layer, 
where the cyclotron energy is equal to the photon energy in 
the rest frame of the electron (e.g., Blandford & Scharlemann 
1975). Two types of models have been discussed for the scat- 
tering region: a stationary nonrelativistic plasma in the closed 
magnetic field lines zone (Rajagopal & Romani 1997; Wang 
et al. 1998; Ruderman 2003) and streams of ultrarelativistic 
electron-positron pairs ejected along the open field lines (e.g., 
Lyubarskii & Petrova 1998, 2000, and references therein). In 
the latter case, the effects of the resonant inverse Compton 
scattering on the properties of observed UV/X-ray radiation 
have not been investigated in detail; however, crude estimates 
show that an extremely large pair multiplicity is needed to 
reach an optical thickness of > 1 . In the case of nonrelativis- 
tic plasma in the closed zone, which can be supported against 
the gravitational force by the thermal radiation pressure en- 
hanced by the cyclotron resonance (Mitrofanov & Pavlov 
1982; Rajagopal & Romani 1997), the effects of resonant cy- 
clotron scattering become significant if the electron/positron 
number density is a factor of ~ 10 2 larger than the corota- 
tion (Goldreich-Julian) density, «gj ~ 10 13 cm -3 forGeminga. 
The electron-positron pairs could be supplied from acceler- 
ation zones, but rapid pair production (large multiplicity) is 
needed to provide so high densities. In addition, it remains 
unclear how the electrons/positrons would lose the longitudi- 
nal momentum to become nonrelativistic particles (the trans- 
verse momentum is essentially nonrelativistic due to the fast 
synchrotron/cyclotron losses). If, nevertheless, there is such 
a nonrelativistic plasma screen in the closed zone, the wave- 
length dependence of its optical thickness depends on spatial 
distribution of scattering particles. In particular, the assump- 
tion that the minima in the UV light curves of Geminga are 
caused by such a rotating screen implies a significant amount 
of electron-positron pairs at a distance of ~ 15/?ns, where 
the magnetic field is <~ 10 9 G. One might speculate that the 
broader minima in the soft X-ray light curve are also caused 
by a partial eclipse by the screen. In this case, the X-ray res- 
onance layer (at a distance ~ 5/?ns) should subtend a larger 
solid angle than the UV resonance layer. Alternatively, if the 
soft X-rays pulsations are caused not by the screen but by 
the intrinsic anisotropy of the thermal radition in the strong 
magnetic field (see above), then the magnetosphere is trans- 
parent for the X-rays we observe, i.e. the inner boundary of 
the plasma screen is located beyond <~ 5/?ns- 

To infer the size of the putative screen and understand the 
spatial distribution of electrons/positrons, independent infor- 
mation on the orientation of the spin and magnetic axes would 
be very useful. Such information could be obtained from X- 
ray pulsations of the TH component emitted by hot polar caps, 
which apparently contributes to the X-ray spectrum around 
~0.7 keV (see Fig. 10). Caraveo et al. (2004b) have inter- 
preted the EPIC data as displaying a varying TH component 
arising from a rotating polar cap. However, in our analy- 



sis the contribution of this component to the spectral flux is 
small and is not clearly seen in the pulse profile. Indeed, 
the 0.6-1.0 keV pulse profile (Fig. 12) can be decomposed 
into a broad thermal pulse with a maximum at <fi « 0.8, sim- 
ilar to the mainly thermal 0.2-0.5 keV pulse, and a narrower 
nonthermal pulse centered at the same phase, </> w 0.6, as the 
higher of the two purely nonthermal peaks in 2-8 keV light 
curve. In other words, we see no convincing evidence for the 
"TH component" (hence, polar caps) in the energy-dependent 
light curves. More definitive information on the axis orienta- 
tion could be obtained from phase-resolved X-ray polarimetry 
(e.g., Pavlov & Zavlin 2000; Lai & Ho 2003), but it is not pos- 
sible with the currently operating X-ray missions. 

4.2. Nonthermal emission 

The multiwavelength observations of Geminga show that its 
emission in the NIR-optical, hard X-rays, and 7-rays is non- 
thermal, presumably generated in the NS magnetosphere. The 
comparative analysis of the results of these observations pro- 
vides an opportunity to understand the mechanisms responsi- 
ble for the magnetospheric emission in different energy bands. 

First, we can compare the phase-integrated X-ray spectra 
at E > 2 keV and optical spectra at E < 3 eV, where the 
nonthermal emission dominates. As shown in §3.1, the fits 
to the X-ray spectrum yield substantially different slopes of 
the PL component, depending on whether the TH component 
is included in the model. With the data available, we can- 
not statistically prove or reject the TH component. However, 
we can extrapolate the PL components of the two different 
models (with and without the TH component) to the opti- 
cal and compare the extrapolations with the observed non- 
thermal spectra (Figs. 7 and 10). The PL fit of the 2.5-10 
keV tail has about the same slope as the optical PL compo- 
nent [T x = 1.56 ±0.24, r = 1.46 ±0.12 and 1.41 ±0.13 for 
E(B-V) = 0.03 and 0.07, respectively], and the extrapolation 
of the X-ray PL spectrum is marginally consistent with the 
optical fluxes (the uncertainty of the extrapolation is shown 
in Fig. 10 [lower panel] and Fig. 14 for absorbed and unab- 
sorbed spectra, respectively). The NIR through X-ray spec- 
trum can be crudely described by a PL model with a slope 
T x ~ 1-3, which is smaller than T x and T , but the dif- 
ferences are close to 1 a uncertainties of the photon indices. 
Similar NIR through X-ray behavior, with Fx ~ To w 1 .5, has 
been seen in another middle-aged pulsar, B0656+14 (Pavlov 
et al. 2002; Zavlin & Pavlov 2004a). Such behavior implies 
that the optical and X-ray emission are generated by the same 
population of relativistic particles with a PL energy spectrum, 
and by the same (likely synchrotron) mechanism. This con- 
clusion is supported by the fact that the ratio of optical-to-X- 
ray luminosities is about the same for all the pulsars observed 
in both X-rays and optical, despite a large scatter of the X-ray 
and optical "efficiencies", L x /E and Ln/E (Zavlin & Pavlov 
2004b). 

If the TH component is not included in the fit, the extrapo- 
lation of the PL component (T x = 2.02 ± 0.05) of the BB+PL 
fit of the 0.2-10 keV spectrum exceeds the optical fluxes by 
more than two orders of magnitude (see the upper panel of 
Fig. 10). This might imply a flattening of the spectrum with 
decreasing energy, as observed in the younger Crab and Vela 
pulsars (Sollerman et al. 2000; Romani et al. 2005). However, 
the continuation toward higher energies of the NIR-optical 
PL component crosses the X-ray PL component at energies 
> 10 keV, which means a "double break" of the spectrum be- 
tween the optical and X-ray bands (i.e., T becomes smaller 
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FIG. 14. — Multiwavelength phase-integrated spectrum of Geminga. The 
EGRET, COMPTEL, and OSSE points are from Mayer-Hasselwander et al. 
(1994), Kuiper et al. (1996), and Strickman et al. (1996), respectively. Fits 
with the PL model for 3 spectral bands (with ±1(7 uncertainties) and their 
extrapolations are shown. 



than T and then increases to I^). Although such behavior 
has been suggested (with a considerable uncertainty) for the 
young LMC pulsar B0540-69 (Serafimovich et al. 2004), it 
was not observed in other pulsars and currently does not look 
very plausible. 

One can also extrapolate the best-fit X-ray PL components 
toward higher energies to compare with the 7-ray data (Fig. 
14). The CGRO EGRET 7-ray spectrum has a slope T 7 = 
1.50 ±0.08 in the energy range 70 MeV to 2 GeV (Mayer- 
Hasselwander et al. 1994), close to that inferred from the PL 
fit of the 2.5-10 keV tail. However, the extrapolation of the 
EGRET spectrum into the X-ray range exceeds the X-ray PL 
tail by about one order of magnitude, while the extrapolation 
of the X-ray spectrum into the EGRET range underpredicts 
the observed flux by a factor of ^20 at 2 GeV. (The discrep- 
ancy is of course much larger if we take the PL component 
of the BB+PL fit.) Interestingly, the extrapolation of the PL 
that crudely connects the optical and X-ray points (Tox ~ 1 -3) 
predicts approximately correct 7-ray fluxes at ~ 0.5 GeV, but 
the difference in the slopes is a factor of 2.5 larger than the 1 a 
error of T 7 . Thus, unless there is a significant systematic error 
in the EGRET data analysis 12 , we have to conclude that there 
should be a double break in the spectrum between the EPIC 
and EGRET bands, which likely means that the 7-ray ray 
emission is generated by a different mechanism (e.g., curva- 
ture radiation), or by a different population of relativistic par- 
ticles, than the optical/X-ray emission. The CGRO COMP- 
TEL and OSSE observations were not sensitive enough to 
prove or reject the existence of such a break (Kuiper et al. 
1996; Strickman et al. 1996). 

In the above discussion on the connection between the NIR- 

12 We should note in this respect that Grenier, Hermsen, & Henriksen 
(1993) reported a substantially softer spectrum from the COS-B observations 
of Geminga: T 7 = 1.84 ± 0.05 in the 0.05-5 GeV range, or 2.02 ± 0.07 in 
0. 14-5 GeV range, with an indication of a spectral turnover below ~ 0.2 GeV. 
With such a soft 7-ray spectrum, a double spectral break in between the EPIC 
and EGRET bands is certainly required. On the other hand, Fierro, Michel- 
son, & Nolan (1998) found a somewhat harder spectrum from the EGRET 
data: r 7 = 1.42 ±0.02 in 0.03-2 Gev range. Given the large uncertainties in 
Fx (including systematic ones: see the discussion in the next paragraph), the 
necessity of the double break becomes less certain. 



optical and hard X-ray phase-integrated spectra, it was as- 
sumed that each of them can be adequately described by a 
PL model. In fact, it is quite plausible that the nonthermal 
high-energy radiation of pulsars is comprised of several com- 
ponents that are peaked in different directions (phases) and 
have different spectral slopes, as it has been recently observed 
in the hard X-ray emission of the Vela pulsar (Harding et al. 
2002) and is routinely seen in the 7-ray range (e.g., Fierro et 
al. 1998). In this case one should expect the phase-integrated 
spectrum to have a concave shape (in a log-log scale), with 
excesses at low and high energies dominated by the softer 
and harder components, respectively. Some indications of 
such a concave spectrum are seen in the Geminga's nonther- 
mal X-ray component: while fitting the high-energy tail of 
the Geminga's EPIC-pn spectrum with a single PL and dif- 
ferent low-energy cut-offs (E m i n ), we noticed that the best- 
fit r increases with decreasing E m { n (e.g., T = 1.45 ±0.30, 
1.56 ±0.24, and 1.60 ±0.20 at £ min = 3.0, 2.5, and 2.0 keV, 
respectively). Although the change of V is not statistically 
significant, such a trend suggests that the tail is concave rather 
than flat, which, in turn, hints that there may be several emis- 
sion components (unless the high-energy tail of the thermal 
component is harder than in the thermal models we used). 
Unfortunately, the S/N at E > 2 keV is too low to verify this 
directly with phase-resolved spectroscopy. 

It is also interesting to compare the nonthermal pulsations 
in different energy bands. The 7-ray light curve (lower panel 
of Fig. 12) shows two peaks per period, at <j> ~ 0.55 (Peak 
1) and <f> ~ 1.05 (Peak 2), with a bridge between them. In 
the X-ray 2-8 keV band, we also see two peaks, but they are 
broader, their separation (A<j> ~ 0.35 or 0.65) differs signif- 
icantly from the A<j) ~ 0.5 in 7-rays, and about 60%-70% of 
the 2-8 keV emission is unpulsed, in contrast to about 100% 
pulsed 7-ray emission. We cannot determine from these data 
alone what is the correspondence (if any) between the X-ray 
and 7-ray peaks (e.g, the higher X-ray peak could correspond 
to either Peak 2 or Peak 1, which would mean that it trails 
the corresponding 7-ray peak by S(f> « 0.55 or 0.05, respec- 
tively). It is tempting to identify the stronger X-ray emis- 
sion at phases ~ 0.1-0.5 as a bridge similar to that in the 
7-ray light curve (which would mean that the higher X-ray 
peak corresponds to Peak 2), but the low S/N at these phases 
makes this identification rather uncertain. Moreover, since no 
statistically significant change of the peak phases with energy 
is seen in the 0.03-3 GeV EGRET range, it is quite possi- 
ble that the 2-8 keV X-ray peaks do not directly correspond 
to the 7-ray peaks, being produced by a distinct mechanism. 
Observations at intermediate energies, 10 keV - 10 MeV, are 
needed to clarify this issue. 

As we mentioned above, the highly asymmetric pulse pro- 
file at intermediate X-ray energies 0.6-1 keV can possibly be 
decomposed into two components: a broad, likely thermal, 
component, centered at <j) « 0.75, similar to that seen in the 
0.2-0.5 keV band, and a narrower nonthermal peak around 
the phase <f) w 0.6 of the higher of two 2-8 keV peaks. How- 
ever, the putative nonthermal peak at 0.6-1 keV looks sig- 
nificantly broader than the corresponding peak at 2-8 keV, 
while the lower 2-8 keV peak is hardly seen in the 0.6-1 keV 
band. The 0.2-0.5 keV light curve is clearly dominated by the 
very broad, thermal peak centered at <f) ~ 0.8; in the simple 
dipole geometry of magnetic field (and the corresponding ax- 
isymmetric temperature distribution), this phase corresponds 
to the nearest approach of the only visible magnetic pole to 
the center of stellar disk. Small contributions from the non- 
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theramal radiation are likely seen even at these low energies: 
e.g., some excess emission over the smooth pulse profile at 
4> ~ 0.55-0.65 might be due to the nonthermal peak centered 
at 4> = 0.6, while the narrow "sub-peak" at (f> w 0.9 might be 
associated with the lower of two 2-8 keV peaks. So, it seems 
that the two nonthermal X-ray peaks persist through the whole 
observed X-ray range, being most clearly seen in the 2-8 keV 
band, where they are not "contaminated" by the thermal radi- 
ation. 

The connection between the optical-UV and X-ray/7-ray 
pulsations (see Figs. 8 and 12) remains unclear. We do see 
two peaks in the FUV band, with about 0.45 phase separa- 
tion, but they may not be related to the hard X-ray and/or 7- 
ray pulsations because the FUV spectrum is mostly thermal. 
We see only one broad peak in the NUV light curve, where 
the nonthermal contribution is substantial, ~ 30%^l-0%, but 
we cannot rule out the possiblity that it is composed of two 
peaks because the NUV light curve is very noisy. Two peaks 
with a separation of about 0.5 in phase, similar to the 7-ray 
peaks, were apparently seen in the optical B band (Shearer et 
al. 1998), but those observations were compromised by a very 
high variable background. 

It is interesting to compare the multiwavelength nonther- 
mal radiation of Geminga with theoretical models. For in- 
stance, the outer gap model by Zhang & Cheng (2001) pre- 
dicts one nonthermal X-ray peak coinciding in phase with the 
7-ray Peak 1 and one soft thermal X-ray peak centered at the 
phase of the 7-ray Peak 2. However, we see two nonther- 
mal X-ray peaks, while the presumably thermal X-ray peak 
is centered in the middle between Peak 1 and Peak 2. The 
latter property is expected in the polar cap models for high- 
energy radiation, which assume that both nonthermal X-rays 
and 7-rays are emitted along a hollow cone inscribed in the 
surface formed by the last open field lines (Harding & Mus- 
limov 1998). In this model, the two nonthermal peaks are seen 
at the phases when the line of sight is tangent to the cone at 
the site where the radiation is generated (a few NS radii from 
the surface), which implies that only one pole is visible (i.e., 
the angles between the rotaion and magnetic axes and rotation 
axis and line of sight are < 30°). The different separation be- 
tween the X-ray peaks, in comparison with that of the 7-ray 
peaks, might mean that they are generated at a different dis- 
tance from the NS (e.g., closer to the NS if the higher of the 
two 2-8 keV peaks corresponds to the 7-ray Peak 1) while 
their larger widths suggest that the X-rays are less beamed 
than the 7-rays. We might speculate that the two peaks be- 
come even broader and closer in the NUV band, so that they 
appear as one peak, but such a trend is not seen in the optical. 

Alternatively, the soft X-ray pulsations may still be mag- 
netospheric in origin. This is, in fact, natural in the outer 
magnetosphere picture of Romani and Yadigaroglu (1995). In 
their model the closest polar cap approach is 5<f) ~ 0.1 before 
the 7-ray Peak 1, while the cap producing the observed 7- 
ray emission passes 6<p = 0.5 later, at (j> w 0.95 in Figure 12. 
These phases are close to the primary and possible secondary 
UV minima of Figure 8. In the Dyks and Rudak (2003) two- 
pole model, the 7-ray peaks arise from opposite hemispheres, 
but the phase of closest pole approach is rather similar. In 
either case, it seems natural to interpret these minima as the 
result of scattering screens above the polar caps removing flux 
(§4.1.2). It would be very surprising if the FUV pulse mini- 
mum level represented the true Rayleigh-Jeans surface flux 
(with the broad peaks being hot cap or non-thermal emission), 
as this would require a remarkably low NS temperature. In 



this interpretation there would be a distribution of soft non- 
thermal components dominating the X-ray pulse, mimicking 
here the TH component and perturbing the fit of the TS com- 
ponent. Little if any of the pulsed emission would come from 
a thermal cap. The (unpulsed) X-ray surface emission would 
be about 30% lower than that fitted for the TS component, al- 
lowing better agreement with the Rayleigh-Jeans UV flux in 
Figure 9. Higher statistics FUV phase-resolved spectra are re- 
quired to determine whether the closest approach of the mag- 
netic axis to the line of sight is associated with the FUV min- 
ima or the soft X-ray maxima. 

5. SUMMARY 

The results of our work can be summarized as follows. 

1. The STIS/MAMA observations of Geminga have al- 
lowed us to detect, for the first time, its FUV radiation, mea- 
sure the FUV spectrum and NUV flux, and detect FUV and 
NUV pulsations. We also measured the flux in the ACS/WFC 
F555W band. To understand the multiwavelength properties 
of the Geminga's radiation, we analyzed its X-ray spectrum 
and pulsations observed with XMM-Newton. 

2. The phase-integrated NIR through FUV spectrum con- 
sists of two components, thermal and nonthermal. The ther- 
mal component, which dominates at A < 3000 A, is emitted 
from the NS surface. Its flux corresponds to a brightness 
temperature T R] = 0.3-0.4 (c/200/^13) 2 MK. The nonthermal 
component, which is likely generated in the pulsar's magneto- 
sphere, shows a PL spectrum with photon index To = 1.3-1.6 
and dominates at NIR-optical wavelengths. We do not con- 
firm the previously claimed spectral feature in the V band. 

3. The phase-integrated X-ray spectrum is dominated by 
a thermal soft component at E < 0.5 keV, with a blackbody 
temperature T s « 0.5 MK and radius R s « 13^200 km. Extrap- 
olation of this spectrum into the optical-UV domain slightly 
overpredicts the thermal component of the observed UV spec- 
trum, in contrast to all the other NSs observed in both X-ray 
and optical-UV ranges. This might be associated with the 
very low temperature or a different chemical composition of 
the Geminga's surface. Under the assumption that the black- 
body model properly describes the thermal spectra while the 
nonthermal spectrum is a simple power law, an additional 
thermal component is required to fit the X-ray spectrum, with 
a higher temperature, 7), « 2 MK, and a small size of the emit- 
ting region, R/, « 50 m. It is not clear whether this compo- 
nent is real or it appears because the simplified model spectra 
were used in the fits. Indeed, several soft non-thermal com- 
ponents could contribute to the X-ray emission, leaving the 
residual thermal surface flux in better agreement with the UV 
Rayleigh-Jeans emission. 

4. The slope of the X-ray nonthermal component, which 
dominates at E > 1 keV, is not well constrained, T x ~ 1 .2-2.0, 
because of a large background at these energies and, possibly, 
some deviations of the spectrum from a simple PL model. Its 
extrapolation into the NIR-optical domain is either marginally 
consistent with the observed fluxes or goes above these fluxes, 
requiring a flattening of the nonthermal spectrum with de- 
creasing photon energy. Most likely, the optical radiation and 
the X-rays are emitted by the same population of relativis- 
tic particles in the Geminga's magnetosphere and generated 
by the same mechanism. On the other side, the continua- 
tions of the nonthermal X-ray and 7-ray spectra do not match 
smoothly, which suggests different radiation mechanisms in 
these energy bands. 

5. Unexpectedly, we detected strong pulsations in the FUV 
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band. The light curve of the predominantly thermal FUV ra- 
diation, with a narrow, deep minimum and possibly another 
minimum shifted by about half a period, differs significantly 
from the light curve of the thermal soft X-ray radiation, which 
shows one broad, smooth peak per period, possibly distorted 
by a small contribution from the magnetospheric component. 
Because the FUV pulsations can hardly be explained by an 
anisotropic temperature and/or magnetic field distributions, 
we suggest that they could be associated with a resonance 
scattering of the thermal UV photons in the NS magneto- 
sphere. 

6. The light curve of the nonthermal X-ray emission shows 
two peaks, as the 7-ray light curve, but the X-ray peaks are 
substantially broader, and their separation differs significantly 
from the half-period separation of the 7-ray peaks. At least 
one of the nonthermal peaks is apparently seen at lower X-ray 
energies, superimposed on the broader peak of thermal emis- 
sion. We see no clear connection between the nonthermal X- 
ray pulsations and single-peaked NUV pulsations. Two peaks 
have been apparently seen in the optical, but this very noisy 
light curve is badly in need of confirmation. 

To conclude, the HST STIS/MAMA observations have 
allowed us to study the UV spectrum and pulsations of 
Geminga. Combined with the results of the previous optical 
and 7-ray observations and our analysis of the XMM-Newton 
data, this study has provided a connection between different 
energy bands for the thermal and magnetospheric components 
and elucidated the multiwavelength picture of the Geminga's 
radiation. However, some important properties of this radi- 
ation remain unclear. First, the true nature of the UV (par- 
ticularly FUV) pulsations is still uncertain. To understand it, 



phase-resolved spectroscopy of the FUV radiation would be 
particularly useful, which would require deeper FU V-MAMA 
observations, quite feasible with the HST STIS if it is brought 
back to life in a future servicing mission. Second, it remains 
unclear whether we indeed see a small, high-temperature po- 
lar cap or this is an artifact caused by the use of simplified 
spectral models. To answer this question we need realis- 
tic models for thermal and nonthermal X-ray radiation and 
phase-resolved spectroscopy at E > 0.5 keV with high S/N. 
Third, we do not understand the connection between the non- 
thermal X-ray emission and 7-ray emission. This problem 
could be resolved by observations in a 10 keV - 10 MeV band 
with a future mission equipped with detectors more sensitive 
in this energy range than the detectors of the RXTE and IN- 
TEGRAL missions. NuSTAR, for example, should be able to 
offer sensitive measurements in the hard X-ray range. 
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